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A recently proposed method based on dispersion theory, that allows to extract the scattering 
length of a hadronic two-body system from corresponding final-state interactions, is generalized to 
the situation where the Coulomb interaction is present. The steps required in a concrete practical 
application are discussed in detail. In addition a thorough examination of the accuracy of the 
proposed method is presented and a comparison is made with results achieved with other methods 
like the Jost-function approach based on the effective-range approximation. Deficiencies of the 
latter method are pointed out. The reliability of the dispersion theory method for extracting also 
- , the effective range is investigated. 
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The scattering length provides not only an important measure for the strength of the interaction in a specific 
hadronic two-body system p| but often allows to draw further more general and thus even more interesting conclusions. 
For example in the case of the proton-proton and neutron-neutron systems the corresponding scattering lengths in 
. the 1 Sq partial wave provide a very sensitive test of charge symmetry in the strong interaction |2|. SU(3) symmetry 
' can be tested by comparing the 1 Sq scattering lengths in the neutron-proton and T, + p systems, which should fulfil the 
C***'" ! relation a np = a s + p in case SU(3) symmetry holds rigorously In the chiral limit the ttN S-wave scattering lengths 
vanish and therefore any deviation from that value is a direct measure for how strongly this symmetry is broken. Here 
especially the isoscalar component is of high interest due to its close link to the sigma term of the nucleon Q . On a 
more phenomenological level the nN scattering length is interesting because the magnitude of its real part is directly 
linked with speculations about the existence of ?7-mesonic hadronic bound states such as n 3 He 0, IE IS IS US E3 
or?7 4 iJe[ll. 

Unfortunately, a direct determination of the scattering length is only feasible in a few cases. It can be done with 
scattering experiments sufficiently close to the reaction threshold so that the effective range expansion can be utilized 
1 , ■ for extracting the scattering length. But in practice such experiments are only possible for charged and also (quasi) 
stable particles, as it is the case with, e.g., pp, n + p or K + p scattering. One can also extract the scattering length 
from a study of the hadronic level shifts of atoms 01 like ^ Pi Tt~d or also pp ^i). F° r the majority of the hadronic 
two-body systems information about the scattering length is only accessible via an investigation of the final-state 
^ . interaction of systems which have at least three particles in the final state. 

While in the former type of experiments the accuracy of the scattering length is directly connected with the 
precision of the data more detailed and often sophisticated considerations are necessary in order to estimate additional 
uncertainies that arise when the scattering length is extracted from final-state effects [l5|]- However, in some cases such 
an estimation is facilitated by the fact that the reaction mechanism is known. For example, the reaction nd — > (nn)p, 
one of the prime sources of the nn scattering length, can be analysed by means of rigorous Faddeev calculations 
[lil IrH ITsf . Reactions involving the pion such as nd, iT~d — ► jnn or 7c? — > n + nn, which can be used to extract the 
ttN and nn scattering lengths, respectively, can be tackled by chiral perturbation theory in a well controlled way in 
the relevant near-threshold regime [T^. l2fjj . 

In a recent publication jgjfl we argued that also large- momentum transfer reactions such as pp — > K + pA [5S US 01 
or jd — > K + nA |5S I2H l27l l28l 123 . l3Cj are excellent candidates for extracting information about the scattering 
lengths. In reactions with large momentum transfer the production process is necessarily of short-ranged nature. As 
a consequence the results are basically insensitive to details of the production mechanism and therefore a reliable 
and general error estimation can be given. Indeed, in |2l| a formalism based on dispersion theory was presented 
that relates spectra from large-momentum transfer reactions, such as pp — > K + pA or 7c? — > K + nA, directly to the 
scattering length of the interaction of the final state particles. The theoretical error of the method was estimated 
to be 0.3 fm or even less, which is comparable to the error quoted in the context of the determination of a nn |3l|. 
say. This estimate was confirmed by comparing results obtained with the proposed formalism to those of microscopic 
model calculations for the specific reaction pp — > K + Ap. The arguments but also the formalism of Ref. |2l| are, of 
course, valid for any production or decay process which is of short-ranged nature, i.e. also for investigation of hadronic 
two- particle subsystems resulting from the decay of the J/Sfr or B mesons [32l l33j . 
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In the present paper we want to investigate further aspects of extracting scattering lengths from final-state inter- 
actions which were not addressed in our earlier work. One of those topics is the presence of the Coulomb interaction. 
In many interesting hadronic two-particle systems both particles carry charges, like in the already mentioned T, + p 
channel whose scattering length could be extracted from the reaction pp — > K°Y, + p. Then the production amplitude 
acquires additional singularities, due to the long-range nature of the Coulomb forces, and the formalism developed in 
Rcf. [2l| is no longer directly applicable. We will derive the modifications that are necessary in order to adapt the 
dispersion-relation method to the situation when the Coulomb force is present in the final-state interaction. We also 
demonstrate in a toy model calculation how one has to proceed in a concrete application to data. 

In addition we present a more detailed examination of the accuracy of the method proposed in Ref. [2jJ . A test 
based on one specific model calculation, namely for the reaction pp — ► K + hp, has been already performed in that 
paper. However, here we want to put this investigation on a broader basis by considering final-state interactions of 
varying strengths, corresponding to a much larger range of values of the scattering length. In addition we take a 
look at the effective range r e as well which can be also extracted by the proposed dispersion-integral method. For 
the effective range a sensible error estimation is not possible, as was already pointed out in Ref. |2J|, but it is still 
interesting to examine in concrete applications whether meaningful results could be achieved. Finally, and equally 
important, we want to compare the present method with the performance of other, approximative treatments of 
the final-state interaction that are commonly used in the literature to extract information on the scattering length 
and also the effective range. This concerns in particular the Jost-function approach 34] based on the effective-range 
approximation and an even simpler approach that relies simply on utilizing the effective range approximation itself 
|35| . Thereby, we will show that the latter methods lead to (partly drastic) systematic deviations from the true values 
and therefore one has to be rather cautious in the interpretation of results achieved with those methods. 

The paper is structured in the following way: In the subsequent section we give a short review of the dispersion 
integral method for extracting the scattering length from final-state interactions. In section 3 results of an examination 
of the accuracy of this method are presented. Thereby, we consider various (singlet and triplet) S wave YN and NN 
interactions and compare the scattering lengths and effective range extracted with the dispersion-integral method from 
appropriately generated final-state effects with the ones predicted by the models. We also apply two approximative 
methods for treating final-state effects, namely the Jost-function approach based on the effective range approximation 
(Jost-ERA) as well as the effective range approximation itself, and compare their performance with the one of our 
method. In section 4 we generalize the dispersion-integral method to the case where a repulsive Coulomb interaction 
is present in the final state. Test calculations for a final-state interaction with Coulomb are then presented in sect. 5 
and it is discussed in detail how one has to procede in a practical application. The paper ends with a short summary. 



II. FORMALISM 



Our method, which goes back to an idea of Geshkenbein 36, 37], is based on using the dispersion relation technique. 
Consider the production amplitude As of a 2 — ► 3 reaction. To be concrete we discuss pp — > K + pA, or jd — > K + nh, 
with the A7V system being in an L = partial wave and a specific spin state S ( 1 S'o or 3 <Si). This amplitude depends 
on the total energy squared s = (pi +P2) 2 , the invariant mass squared of the outgoing AN system m 2 = (pn+Pa) 2 and 
the momentum transfer t — (p\ — Pk+) 2 , where px, P2, Pn, Pa, and p K + are the 4-momenta of the two initial particles, 
final nucleon, lambda, and kaon, respectively. Then one can write down a dispersion relation for this amplitude with 
respect to m 2 at fixed s and t 

A s (s 1 t,m) = - —dm + - I ^— dm , (1) 



where m 2 is the upper boundary of the lefthand cut, m\ = (mjv + toa) 2 , and 

D s {s, t, m 2 ) = -(A s (s, t, m 2 + iO) - A s (s, t, m 2 - »0)) (2) 

is the discontinuity of the amplitude along the cuts. We neglect here the contributions from possible kaon-baryon 
interactions. In case they are not small, they still can be considered as constant (weakly mass dependent) if one 
chooses the kinematics such that the excess energy of the reaction is significantly larger than the typical range of the 
AN interaction, cf. the discussion in Ref. [2lj . The index S denoting the spin state will be suppressed in the following 
to simplify the notation. 

For a purely elastic A7V system, the discontinuity along the righthand cut would be given by 



D(s,t,m 2 ) = A(s,t, m 2 )e~ iS sinS, 



(3) 
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where 5 is the AN ( 1 5'o or 3 Si) scattering phase shift. Then the solution of Eq. £Q) in the physical region reads (see 

Refs. mmm) 



A(s, t, m 2 ) — exp 



5{m' 2 ) 



2 m' 2 — m 2 — iO 



dm' 



Hs,t,m 2 ), 



(4) 



where $(s,t,m 2 ) contains only lefthand singularities and therefore is a slowly varying function of m 2 . In order to 
ensure this requirement to be fulfilled it is important that the momentum transfer t is large. We assume also that 
there is no bound state in the AN system. 

Consider now a realistic situation where inelastic channels are present - as it is the case with A./V due to the coupling 
to the EiV channel, say. Then one can write down a formula similar to Eq. Q), but with the integration performed 
over a finite range of masses [2lJ : 



A(m 2 ) = exp 



8{m' 2 ) 



ml 2 — m? — iO 



-dm' 



<l(m 2 ) , 



(5) 



where $(m 2 ) is again a slowly varying function of m 2 given the phase shift S is sufficiently small in the vicinity of 



??? 



[21|. The upper limit m ma x has to be chosen in such a way that the corresponding relative momentum of the 
A7V system, p max , is of the order of the typical scale of the A7V interaction, i.e. in the order of 1/a or 1/r. The 
equation 10 can be solved with respect to the S |2l| : 



5{m 2 ) 



2tt 



log \A{m' 2 )/<S>{m 2 max , m' 2 )| 2 / to 2 
\J m' 2 — TOq (m! 2 — m 2 ) 



,/2 



dm' 



(6) 



If one neglects the mass dependence of $(m 2 ) and uses the relation between the partial cross section as and the 
amplitude 



fas 
dm' 2 dt 



(x P '\A s (s,t,m' 2 )\ 2 , 



then one obtains the expression for the scattering length in terms of observables 

. / ' m\ + mN ' 
as — lim — 



m 2 — >m 



1 



\/ m 1 2 — ttIq (m! 2 — m 2 ) 




(7) 



and analogously for the effective range r e . 



III. ACCURACY OF THE METHOD AND COMPARISON WITH OTHER APPROACHES 

The most important advantage of the method proposed by us [2l| is that a reliable estimate for the uncertainty 
of the extracted scattering length can be given. The sources for the uncertainty are trifold: (i) A possible influence 
of the final-state interaction in the other outgoing channels. For the reaction pp — * K + Ap considered in Ref. |2l| 
this concerns the KA and KN systems, (ii) The adopted value for m 2 max , the upper limit chosen for the dispersion 
integral in Eq. (|7|. (iii) A sensitivity to left-hand cuts of the production operation. A detailed analysis of the issues 
(ii) and (iii) , based on general arguments, presented in Ref. [2l| suggests that the error in the scattering length 
should be typically in the order of 0.3 fm or less. The role of issue (i) cannot be quantified theoretically but has to 
be investigated by performing experiments and corresponding analyses at different beam momenta plj . 

In this section we want to present a thorough examination of the accuracy of the proposed method and, in particular, 
to corroborate the error estimate, by means of concrete model calculations. A test based on one specific model 
calculation, namely for the reaction pp — > K + Ap, has been already performed in Ref. |2l|. However, here we want 
to put this investigation on a broader basis by considering final-state interactions of varying strengths, corresponding 
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to a much larger range of values of the scattering length. Furthermore, and equally important, we want to compare 
the present method with the performance of other, approximative treatments of the final-state interaction that are 
commonly used in the literature to extract information on the scattering length and the effective range and have been 
applied to pp — > K + Ap (IHI^. 

One of those approximative treatments follows from the assumption that the phase shifts are given by the first two 
terms in the effective range expansion, 

p C ot(S(m 2 )) = -- + ^p 2 , (8) 
a 2 

usually called the effective range approximation (ERA), over the whole energy range. Here p is the relative momentum 
of the final-state particles under consideration in their center of mass system, corresponding to the invariant mass 
to 2 . In this case the relevant integrals (QJ can be evaluated in closed form as [34( 

,/ 2n (p 2 + a 2 )r e /2 

A(m 2 ) cx - , ; — ^4 , 9 

1 ' -l/a+(r e /2)pi -ip V ' 



where a = l/r e (l + yl — 2r e /a). Because of its simplicity Eq. © is often used for the treatment of the final-state 
interaction (FSI). 

A further simplification can be made if one assumes that a S> r e . This situation is practically realized in the 1 Sq 
partial wave of the NN system. Then the energy dependence of the quantity in Eq. @ is given by the energy 
dependence of the elastic amplitude 

A(m 2 ) oc ; , * n , (10) 

1 ' -l/a+(r e /2)p 2 -ip V 1 

as long as p <C 1/V e . Therefore one expects that, at least for small kinetic energies, NN elastic scattering and meson 
production in NN collisions with a NN final state exhibit the same energy dependence [3 HH > which indeed was 
experimentally confirmed. This treatment of FSI effects is often referred to as Migdal- Watson (MW) approach 35] ■ 
In order to examine the reliability of the three methods described above we took different YN models from the 
literature |43l E3. EH l4fil] and calculated the production amplitude A{m 2 ) utilizing the meson exchange model from 
Ref. |Z3 . Then this amplitude was used for extracting the scattering length by means of the dispersion integral Eq. 
|J7|| or from the approximative prescriptions given by Eqs. JJjJ and IjlQI) . For comparison we considered also the 1 Sq 
partial wave of the np system of the Argonne potential [48| . In this case A(m 2 ) was set equal to the scattering wave 
function ^(p, r) at the origin, more precisely to ty~(p, 0)*, which corresponds to the assumption that the production 
operator is point-like. 



TABLE I: S'-wave scattering lengths a (in fm) for various YN |44l 145^1 and NN |4^l potentials. The results for the original 
models are compared with those obtained by applying the dispersion integral method (7) and the approximations Eq. © 
(Jost-ERA) and Eq. JTDJ (MW). 



model 


exact result 


disp.int. 


Jost-ERA 


MW 


Jiilich 01 singlet 


-1.02 


-1.03 


-1.28 


-1.67 


Nijmegen 97a singlet 


-0.73 


-0.75 


-0.98 


-1.33 


Nijmegen 97f singlet 


-2.59 


-2.57 


-2.96 


-3.35 


Jiilich 01 triplet 


-1.89 


-1.66 


-2.05 


-2.42 


Nijmegen 97a triplet 


-2.13 


-1.98 


-2.37 


-2.75 


Nijmegen 97f triplet 


-1.69 


-1.61 


-2.00 


-2.37 


Argonne vl4 singlet 


-23.71 


-23.54 


-24.56 


-24.79 



Some selective results (for the Nijmegen NSC97 Q and Jiilich 01 YN models and the Argonne vl4 Q NN 
potential) are summarized in Table [I] The second column contains the correct scattering length evaluated directly 
from the potential model. One can see that the extraction of the scattering length via the dispersion integral Q 
yields results pretty close to the original values for all considered potentials. In fact, in most cases the deviation 
is significantly smaller than the uncertainty of the method, estimated in Ref. [2l| to be 0.3 fm. The results of the 
Jost-ERA approach Eq. JjJJ exhibit a systematic offset in the order of 0.3 fm. The situation is much worse for the 
Migdal- Watson approach Eq. (|10|l where a similar offset is found though now in the order of 0.6 fm. As a consequence, 
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the extracted values differ by 50 % or more from the correct scattering lengths. Only for the 1 S'o np partial wave 
the disagreement is still in the order of 5 %. Here the reliability of the Jost-ERA and Migdal- Watson approaches are 
comparable. This is in agreement with the expectations mentioned above. 

The systematic offset inherent in the Jost-ERA approach as well as in the Migdal- Watson prescription can be best 
seen in Fig.^ where we shown the difference between the scattering lengths predicted by various models and the values 
extracted via the dispersion integral (circles), the Jost-ERA method (squares) and the Migdal- Watson prescription 
(triangles). 

While the Jost-ERA approach might still be a reasonable tool for getting a first rough estimate of the scattering 
length for a particular two-body interaction one should be rather cautious when using it for more quantitative analyses. 
In particular, its application in a combined fit to elastic scattering data and invariant mass spectra, e.g. to Ap and 
pp — > K + Ap, is rather problematic and can easily cause misleading results. Because of the offset in the scattering 
length in applications to final-state effects it is clear that a combined fit cannot converge to a unique (the "true" ) Ap 
scattering length. Only the elastic data will favour values close to the "true" scattering length whereas the production 
data tend to support larger (negative) values. This is obvious from the corresponding Jost-ERA results presented 
in Table [I] and also from Fig. ^ We believe that the analysis of Hinterberger and Sibirtsev presented in Ref. ^2] 
is an instructive exemplification of this dilemma. Employing the Jost-ERA approach to low energy total Ap cross 
sections ^|l,|5(j and to experimental results for the missing mass spectrum of the reaction pp — > K + X [ElJ separately, 
they derived (spin averaged) scattering lengths of a = — 1.81 jlo'21 ^ m an( ^ a = ^^-^-0 23 respectively. Taking 
into account the error bars this is roughly the difference we would expect from the offset (of around 0.3 fm) seen 
in our test calculations and, therefore, one must consider the results as being practically consistent with each other. 
But the authors of Ref. ^2 attempted to "reconcile" the results even more by introducing a spin-dependence in the 
fitting procedure. Indeed, with the relative magnitude of singlet to triplet contribution in the production reaction as 
free parameter (their relative strength in the elastic channel is fixed at 1:3 by the spin weight!) a "fully consistent" 
description of the combined data could be achieved and apparently the spin-singlet as well as spin-triplet Ap 
S-w&ve scattering lengths could be determined from spin-averaged observables. Our experience with the Jost-ERA 
approach reported above, however, strongly suggests that the sensitivity to the spin seen in this analysis is most likely 
just an artifact of the method applied. 

TABLE II: S'-wave effective ranges r e (in fm) for various YN |44T.|45| and NN 48] potentials. The results for the original models 
are compared with those obtained by applying the dispersion integral method J7J and the approximation Eq. © (Jost-ERA) . 



model 


exact result 


disp.int. 


Jost-ERA 


Jiilich 01 singlet 


4.49 


4.31 


2.48 


Nijmegen 97a singlet 


6.01 


4.78 


2.81 


Nijmegen 97f singlet 


3.05 


2.82 


1.60 


Jiilich 01 triplet 


2.57 


2.49 


1.89 


Nijmegen 97a triplet 


2.74 


2.60 


1.70 


Nijmegen 97f triplet 


3.34 


2.67 


1.69 


Argonne vl4 singlet 


2.78 


2.91 


0.43 



Let us now come to the effective range r e . Since the dispersion relations yield only an integral representation for 
the product a 2 ((2/3)a — r e ) but not for the effective range r e alone [2l| it follows that the attainable accuracy of r e is 
always limited roughly by twice the relative error on a. Still, it is interesting to see what values one gets for r e from the 
dispersion integrals. Corresponding results are presented in Table Inland in Fig. [3 Evidently, the values extracted via 
the dispersion integral agree much better with the original results as one might have expected. In fact, in practically 
all cases the deviation is in the order of only 5 % or even less. This suggests that one could use the dispersion integrals 
also to extract the effective range r e from data. But one should keep in mind that, unlike the case of the scattering 
length, now one cannot rely on a solid and general estimate of the uncertainty. As far as the Jost-ERA approach is 
concerned it is clear from Table ITT1 that it yields rather poor results. In case of the Migdal- Watson prescription (|10|) it 
turned out that the fit always prefers an effective range r e equal to zero. This is due to the term proportional to r^p 4 
in the denominator of the A(m 2 ) that make the production cross section decrease too fast as compared to the data 
(or to our calculations with realistic models) . Therefore we don't show any results of the Migdal- Watson fit for r e . 

One should note here that the upper limit in the dispersion integrals was always taken such that p m ax = 205 MeV/c 
(as in |21|) that corresponds to e max = m max — m as 40 MeV for the AN and T,N systems and e max « 45 MeV for 
the NN system. In the latter case it is interesting to see what happens if one varies the range of integration, since 
the energy structure in the NN interaction is much narrower due to the large NN scattering length. For e max = 10 
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MeV and e max = 20 MeV as upper limits of the integration one gets the scattering lengths a — —22.62 fm and 
a = —23.17 fm, respectively - which are in principle still close to the original value. For the effective range, however, 
the calculation yields r e = 4.78 fm and 3.72 fm, respectively. The reason why the agreement for e max > 40 MeV is 
so good is that the NN phase shift becomes sufficiently small at such energies, which implies a small uncertainty 
according to the error estimation in Ref . [2lJ . 



IV. DISPERSION RELATION IN THE PRESENCE OF COULOMB REPULSION 

In the case when both baryons in the final state carry charges (for example in the reaction pp — > K°pT, + ) there 
is a Coulomb interaction between them. Then the production amplitude A(m 2 ) acquires additional singularities at 
p = 0, due to the long-range nature of the Coulomb forces, and the formalism developed in Sect. 2 is no longer 
applicable directly. In this section we want to describe the modifications that are necessary in order to adapt the 
dispersion-relation method to the situation when the Coulomb force is present in the final-state interaction. We 
restrict ourselves to the case of a repulsive Coulomb interaction so that no bound states are present. 

In order to elucidate the principle idea we start out from the case of elastic (two-body) scattering. Here the 
problem can be most conveniently dealt with by applying the Gell-Mann-Goldberger two-potential formalism [H^ . 
Let us assume that the total potential V = V c + V S is given by the sum of a short-ranged hadronic potential V s and 
the Coulomb interaction V c . Then the total reaction amplitude T can be written as T = T c + T CSl where T c is the 
Coulomb amplitude and T cs is defined by 



T cs = (l + T c G )t cs (l + G T c ), (11) 

where t cs fulfils a Lippmann-Schwinger equation, 

t cs = V s + V s G c t cs , (12) 

with the short-range potential V s as driving term. To obtain the physical on-shell amplitudes one needs to project the 
corresponding T-operators on the so-called Coulombian asymptotic states |poo±) which are related to the Coulomb 
scattering states (with fixed angular momentum - in our case I — 0) |jJ±) c via |p±) c = \poo±) + GqT^ 1 |poo±) H3- 
Here p denotes the center of mass momentum in the baryon-baryon system. In this way one obtains in particular 

c { v -\t cs \p+) c =-^ 1 f cs {p), (13) 

where f cs is the so-called Coulomb-modified nuclear scattering amplitude and (i is the reduced mass. Its relation to 
the phase shift 5 CS is the following 

I- = % ■ (14) 

2 

with S c denoting the pure Coulomb S-w&ve phase shift given by 5 C — arg (r(l + irf)) with 77 = 

It has been shown in Ref. |54| under rather general assumptions that the modified amplitude f(p) — 
e~ 2%s ° f cs (p) I 'C 2 (p) is free of the Coulomb sing ularities on the physical sheet and possesses only the singularities 
caused by dynamical cuts (see also Refs. |55l l56l l57l IH^ 1 ). In addition below the inelastic cuts and above the 
two-baryon threshold the modified unitarity relation reads 

f(s + iO) - f(s - tO) = 2ipJ{p)~f*{p)C 2 (p) (15) 

with C 2 (p) — -r^zi being the Coulomb penetration factor. 

Furthermore an effective range function, modified for the presence of the Coulomb interaction, can be defined as 
well. Is is given by 

S(p) = pC 2 {p) cot 6 CS ( P ) + Q(p) = -l/o M + r eP 2 /2 + .., (16) 

where Q(p) = [ie 2 [ip{ir]) + ip{-irj) - 2lnrj], ip{z) = T'(z)/T(z). 

Coming back now to the production reaction it can be analogously shown that also the modified production 
amplitude 

A(m 2 ) = e- lS "A(m 2 )/C(p) (17) 
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is free of the aforementioned singularities . Therefore a dispersion relation similar to Eq. JQ) can be written down 

7, 2, 1 r 2 D(m' 2 ) , , 2 1 /- 00 £>(m' 2 ) , , 2 

A{m 2 ) = - —\ %dm' 2 + - / — ^ %dm' 2 . (18) 



Unitarity implies that the discontinuity for the elastic cut is 

D(m 2 ) = A{m 2 )e- lS "° sin 5 CS . (19) 
The solution to Eq. Ijl8[l is found in complete analogy to the case without the presence of the Coulomb interaction, 

'l S cs (m> 2 ) 



A(m ) = exp 



7T 



■o 



i' 2 — m 2 — iO 



dm 



*(m") , (20) 



where ^>(m 2 ) is some function slowly varying with to 2 . If one neglects the the weak to 2 dependence present in ^(to 2 ) 
the expression for the phase shift <5 CS in terms of the differential cross section becomes 



8 cs (m 2 ) _ 1 



log 



1 d 2 a 



p'C 2 (p') dm' 2 dt 



%dm 12 . (21) 



a/to 2 - ml 27r J m 2 ^/to' 2 - TOq (to' 2 - m 2 ) y rn 2 lax - 
Using the effective range expansion (| 1 fcil) one can then extract the scattering length a cs from this dispersion integral. 

V. TEST OF THE METHOD FOR THE COULOMB CASE 

One of the obvious reactions for applying the formalism with Coulomb is pp — » K°Yi + p where one could determine 
the T,N scattering length for the isospin 3/2 state. Note that the S + p channel does not couple to the AiV system 
and is therefore free of inelastic cuts (that start already on the left-hand side) as required for the applicability of 
the disperson integral method. For this reaction one could perform a model calculation analogous to the one for 
pp — > K + Ap 01 which we used for testing the dispersion-integral method in the absence of the Coulomb interaction 
|21j. However, the implementation of Coulomb effects into our momentum-space code is technically complicated and 
requires also some approximations [EjJ. Thus, for the present test calculation we adopt a different strategy. First, 
instead of the momentum-space YN models of Refs. 0, HE H^ we take the r-space Argonne (NN) potential, 
however, with parameters modified in such a way that the effective range parameters are similar to those predicted 
by realistic YN potentials |43l EH liH li^ for the T,N I = 3/2 1 Sq partial wave. In particular we prepared two models 
with Coulomb modified scattering length of a cs — —3.24 fm (model 1) and a cs = —1.86 fm (model 2), respectively. 
The corresponding scattering lengths without Coulomb interaction are —4.11 fm and —2.01 fm, respectively. For the 
transition amplitude we use the scattering wave function calculated from those potential models and evaluated at the 
origin. This corresponds to the assumption that the production operator is point-like, which is reasonable as long 
as we are interested only in the mass dependence of the production amplitude. The corrections stemming from a 
possible mass dependence of the production operator were discussed in Ref. |2ll . 

The results of applying Eq. J7J) with m max — m = 40 MeV (p ma x — 205 MeV/c) are shown in Fig. where 
we plot the function 1/S(p) which should coincide with the scattering length — a cs at p = 0. Obviously, there is 
a strongly nonanalytic behavior of the extracted inverse effective range function when approaching the threshold 
- which, however, can be easily understood. It is clear from Eq. that the threshold behavior of the Coulomb 
modified phase shift is S cs w — a cs pC 2 (p), i.e. S cs goes to zero very rapidly. To obtain such a behavior on the left-hand 
side of Eq. (|21|l one needs to have a very precise cancellation in the integral on the right-hand side of Eq. I|21|l that is, 
of course, impossible if one truncates the integral at a finite momentum. But still one can expect for Eq. I|21|l to work 
for momenta not too close to the threshold, namely above the typical Coulomb scale of 2wa/ ' \i « 25 MeV/c where 
the factors C 2 (p) and Q(p) that appear in the effective range function (Eq. (|16fl ) become smoother. This is indeed 
the case, as can be seen from Fig. [3] Thus, a natural and practical step here would be to extrapolate the extracted 
S(p) to the threshold from above. Using a 4-th order polynomial of the type — l/a cs + r e p 2 /2 — Pr^p 4 and fixing 
the coefficients in the region 50 — 100 MeV/c, i.e. well above the Coulomb structure, one can then extrapolate S(p) 
to the threshold, cf. the dash-dotted lines in Fig. In this case one gets a satisfactory agreement between the true 
and extracted scattering lengths. In fact, the deviations are not worse then in the case when we consider the same 
potentials without Coulomb interaction and they are also within the theoretical error of 0.3 fm estimated in Ref. [2l| . 
The extracted values are a cs = —3.10 fm and a cs = —1.86 fm for models 1 and 2, respectively, with Coulomb, and 
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a s = —4.05 fm and a s — —2.06 fm when the Coulomb interaction is switched off. We also checked that the sensitivity 
of the result to the region of interpolation of the effecitve range function S(p) is rather low. For instance if one shifts 
the lower bound of this region to 70 MeV/c the corresponding change in the scattering length will be less then 0.05 
fm. 

If one wishes to consider a more realistic situation, one needs to deal with mass distributions with finite statistical 
errors, finite mass resolution, and finite binning. To examine also this situation we have generated two data sets, 
corresponding to the models 1 and 2 as shown on Fig. 01 We have chosen the binning as well as the mass resolution 
to be equal to 2 MeV (the same as in the experiment 51] that was analyzed in 21]), and the statistics to be rather 
high in order to minimize the influence of the statistical error bars on the results. The excess energy was set to 40 
MeV to simplify the simulation. In a realistic situation larger values are preferable to minimize the influence of the 
meson-baryon interactions, cf. the corresponding remarks in Ref. |2l| . In our test calculation such meson-baryon 
interactions are neglected anyway. 

We start here with the procedure suggested in appendix A of |2l|, namely by fitting the cross section with an 
exponential parameterization of the type 



d 2 a 
dm 2 dt 



= C 2 (p) x exp 



Co 



G\ 



(m 2 



x phase space . 



(22) 



This formula fits the generated cross section with the \ 2 P er degree of freedom of Xdof ~ 1) c f- Fig. El A new 
problem that arises here is that the production amplitude contains a very narrow structure (of the size 2ira/ 'fi) close 
to threshold as can be seen in Fig. |S] (solid lines). Clearly, this structure can not be reproduced after the fitting 
procedure as it gets smeared out by the mass resolution and binning (their size is much larger than the scale of the 
structure). The amplitudes coming from the fit are depicted in Fig. [S] by dash-dotted lines. However the fit can be 
improved if one notes that the structure comes mostly from the part of the dispersion integral (|20() containing the 
leading term in the S cs expansion near the threshold, namely (in the nonrelativistic case) 



exp 



o P 



-a cs C 2 (p>) fp 2 



P 



l 2 



p^ 



i0 



dp 1 



exp | 



sQ(p)} 



(23) 



where we made a subtraction at p — to render the integral convergent. This does not change the energy dependence 
of the resulting exponent. Indeed the structure disappears after dividing the production amplitude by the factor 
exp[-a cs Q(p)]. 

Obviously the scattering length is unknown before its extraction! But what can be done is to resort to an iterative 
procedure including first the extraction of the unimproved scattering length, then putting it into the fit function, 



dm 2 dt 



C 2 {p) x exp 



Co 



(m 2 



CI) 



2a cs Q(p) 



x phase space 



(24) 



and repeating this step until the procedure converges. Fortunately the iterations converge already after three or four 
iterations, and the resulting amplitudes are shown in Fig. [S] by the dashed lines. The improvement of the fit is quite 
obvious. Finally we applied the combination of the extrapolation and iteration procedures to obtain the scattering 
lengths from the pseudodata. Since the data have a statistical uncertainty we generated a sample of 1000 mass 
distributions and looked at the average value of the scattering lengths. They turned out to be —2.89 ± 0.06 fm for 
model 1 and — 1.82 ± 0.05 for model 2. Note that the deviation from the correct values is now a bit larger, but still 
reasonable ( 0.35 fm in the worst case). Fig. shows how the extracted inverse effective range function approaches 
to the correct one given by the model by the example of model 1. 



VI. SUMMARY 

In a recent publication |2l| we have presented a formalism based on dispersion relations that allows one to relate 
spectra from large-momentum transfer reactions, such as pp — > K + pA or 7c? — > K + nA, directly to the scattering 
length of the interaction of the final-state particles. An estimation of the systematic uncertainties of that method, 
relying on general arguments, led to the conclusion that the theoretical error in the extracted scattering length should 
be less than 0.3 fm. This finding was corroborated in an application of the method to results of a microscopic model 
calculation for pp — ► K + pA. 

In the present paper this dispersion theoretical method was generalized to the case where a repulsive Coulomb force 
is present in the final-state interaction. As an example let us mention the reaction pp — > K pT, + which could be used 
to extract the pT, + scattering length. Though the generalization of the formalism itself is straight forward it turned 



9 



out that there are some additional features due to the Coulomb interaction which need to be taken into account in 
concrete applications of the method to data. These practical aspects were thoroughly discussed and it was shown 
how to circumvent the difficulties. In a test calculation utilizing potential models with effective range parameters 
similar to those of realistic YN interactions the extracted values for the scattering lengths were found to agree within 
0.3 fm with those predicted by the models. Thus, the accuracy of the dispersion theoretical method for extracting 
the scattering lengths from final-state interactions with Coulomb force is comparable to the case where no Coulomb 
interaction is present. 

We presented also a more detailed examination of the accuracy of the dispersion-integral method than in Ref. 
( 12 1 1 ] . In particular we considered final-state interactions of varying strengths, corresponding to a much larger range of 
values of the scattering length. These investigations confirmed the reliability of the general error estimate provided in 
Ref. |2lJ . Indeed in most of the considered cases the deviation of the extracted scattering length from the true value 
was significantly smaller than the uncertainty of 0.3 fm derived in that paper. In addition we studied the effective 
range r e which can be also extracted by the proposed dispersion-integral method. For most of the interaction models 
considered the extracted values of r e agreed remarkably good with the true results. Thus, it might be sensible to use 
the proposed method to extract the effective range from data - though one should always keep in mind that for this 
quantity a generally valid error estimation is not possible |21| . 

Finally, we compared the present method with the performance of other, approximative treatments of the final- 
state interaction that are commonly used in the literature to extract information on the scattering length and also the 
effective range. In particular we tested the Jost-function approach based on the effective-range approximation (Jost- 
ERA) and an even simpler approach that relies simply on utilizing the effective range approximation itself. Thereby, 
we showed that the latter methods lead to systematic deviations from the true values of the scattering lengths in the 
order of 0.3 fm (Jost-ERA) and even 0.7 fm (direct effective-range approximation). This suggests that one should be 
rather cautious in the interpretation of results achieved with those methods. 
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FIG. 1: (Color online) Comparison of different extraction methods for the scattering length a. Shown are the differences 
between results predicted by various YN and NN models and corresponding values extracted via the dispersion integral 
method (circles), the Jost-ERA approach © (squares) and the Migdal- Watson prescription I1U1 (triangles). The lines are 
drawn to guide the eye. 
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FIG. 2: (Color online) Comparison of different extraction methods for the effective range r e . Shown are the differences between 
results predicted by various YN and NN models and corresponding values extracted via the dispersion integral method (circles) 
and the Jost-ERA approach @ (squares). The lines are drawn to guide the eye. 
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FIG. 3: The inverse of the effective range function S(p) for model 1 (upper curves) and model 2 (lower curves). The solid lines 
denote the phase shifts predicted by the corresponding model, whereas the dashed lines correspond to the phase shift extracted 
via Eq. 12 U . The dash-dotted lines show the result of the smooth extrapolation of the dashed lines, as explained in the text. 
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FIG. 4: Pseudo data for the differential E + p cross section generated from model 1 (circles) and model 2 (squares) as a function 
of the Ep invariant mass Me p with corresponding fit by the exponential parameterization Eq. 1221 . 
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FIG. 5: The production amplitude A(rn 2 ) divided by the factor C(p) for model 1 (upper curves) and model 2 (lower curves). 
The solid lines are the amplitudes as calculated from the models. The dash-dotted lines correspond to the fitted amplitudes. 
The dashed lines denote the fitted amplitudes improved by the iterative procedure as discussed in the text. 
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FIG. 6: The inverse of the effective range function S(p) calculated by means of the iterative procedure as discussed in the text. 
The dotted line corresponds to the zeroth iteration, the dashed line corresponds to the first iteration, and the dash-dotted line 
corresponds to the second iteration. The solid curve denotes the exact result. Shown are results for the model 1. 



